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Abstract 

We derive the scale dependence of twist-3 quark-gluon operators, or ETQS matrix elements, at 
one-loop. These operators are used to factorize transverse single spin asymmetries, which are studied 
intensively both in experiment and theory. The scale dependence of two special cases are particularly 
interesting. One is of soft-gluon-pole matrix elements, another is of soft-quark-pole matrix elements. 
From our results the evolutions in the two cases can be obtained. A comparison with existing results 
of soft-gluon-pole matrix elements is made. 



J> , In high energy scattering with a transversely polarized hadron transverse Single Spin Asymme- 

try(SSA) can appear. Such an asymmetry has been observed in various experiments. Reviews about 
this research field can be found in [lj. The study of SSA plays an important role in exploring the inner- 
structure of hadrons. In the framework of collinear factorization, SSA can be factorized with matrix 
| elements defined with twist-3 QCD operators, called ETQS matrix elements [21 [3]. These matrix elements 
contain important information about correlations between more than two partons in the transversely 
polarized hadron and can be extracted from experimental data of SSA. 

To precisely extract the twist-3 matrix elements or predict SSA, one needs to know the scale depen- 
dence of the matrix elements. This dependence plays a similar role like that of standard twist-2 parton 
distributions functions. The later is given by the famous DGLAP equations. In this letter we study the 
scale dependence of these twist-3 matrix elements. Especially, we focus on the matrix elements defined 
with twist-3 quark-gluon operators. 

The definitions of the interested twist-3 matrix elements can conveniently be given with a light- 
cone coordinate system, in which a vector is expressed as a M = (a + ,a~,a±) = ((a + a 3 )/\/2, (a — 
a 3 )/\/2, a 1 , a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . In the system we introduce two light-cone vectors: = (0,1,0,0) 
and V 1 = (1,0,0,0). The anti-symmetric tensor in the transverse space is defined as = e a ^^ v l a np 
with the convention e 0123 = 1. There are 6 twist-3 operators or 6 twist-3 matrix elements for a transversely 
polarized hadron with the momentum = (P + , 0,0,0) and the spin vector = (0, 0, Two of 

them are defined with product of two quark operators and one gluon field strength operator. They are 
defined as [2 [3]: 

2>(si,s 2 ,m) = StfsJ ^^e- lP+ ^^-^) + ^ :El )(P, S |^(y 1 n)7 + G + ^(y 2 n)^(0)|P, S ), 
T A , F ( Xl ,x 2 ) = i s ^|^^ e -^ + (^-*i)+^ (1) 
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with n = 1,2 and s M = e^ v s u . The above definitions are given in the light-cone gauge n ■ G = 0. In other 
gauges, gauge links along the direction n should be supplemented to make them gauge invariant. From 
general principles one can derive the following properties: 

T F (xi,x 2 ) = T F (x 2 ,x 1 ), Ta,f(xi,x 2 ) = -Ta,f(x2,xi). (2) 

The remaining four twist-3 matrix elements are defined only with gluon field strength operators. 
These purely gluonic matrix elements have been given in [5] . Two of the four matrix elements are defined 
with the £[/(iVc')- gauge group constant f abc : 

■{P,s\G a ' +a ( yi n)G b ^(y 2 n)G c ^mP,s), 



T£> G (xi,x 2 ) = -s^ 



(y2(x2-xi)+yixi) 

P+ J 4vr " 
■(P,s\G a ' +a ( yi n)G b ^(y 2 n)G c ^(0)\P, S ). (3) 



The definitions of other two T G and T A G are obtained by replacing if abc with d abc . Again, the above 
definitions are given in the n ■ G = gauge. For the matrix elements with f abc one has: 

T G f \x Xl x 2 ) = -T G f \-x 2 ,- Xl ), T G f \ Xl ,x 2 )=T G f \x 2 , Xl ), 

T^ G ( Xl ,x 2 ) = T^ G (-x 2 ,- Xl ), T { £ G (x 1 ,x 2 ) = -T i £ G (x 2 ,x 1 ). (4) 

Similar relations for T^P and T^ G can also be derived. The defined six twist-3 matrix elements depend 
on the renormailzation scale [i. 

The support of the defined matrix elements can be analyzed with translational covariance. One easily 
obtains that the matrix elements are only nonzero for |a?i 2I < 1 an d \x\ — x 2 \ < 1. There are two 
special cases which are particularly interesting. One is with x\ = x 2 . The corresponding nonzero matrix 
elements are called soft-gluon-pole matrix elements, which describe the correlation between partons and 
one gluon with zero momentum fraction. Another case is with x\ = or x 2 = 0. The corresponding 
matrix elements are called soft-quark-pole matrix elements. The matrix elements in two cases are of 
particular interesting[3l Hj. The scale dependence of soft-gluon-pole matrix elements has been studied in 
[9] . But the obtained results are not completely in agreement. The renormailization of the above 
twist-3 operators in space-time has been studied in [10]. The evolution of twist-3 quark-gluon operators 
has been studied in [11] with the emphasis on the solutions of evolution equations. In this letter we study 
the scale dependence of twist-3 quark-gluon matrix element Tp(xi,x 2 ) and Ta,f{x\, x 2 ) with x± t2 > 
with a different method. From our results one can obtain the scale dependence in the special case of 
x\ = x 2 or of x\ j2 = 0. 

To study the scale dependence of the matrix elements defined with twist-3 quark-gluon operators, it 
is convenient to introduce as in [7]: 

T±(xi,x 2 ) = T F (xi,x 2 ) ± Ta,f{x\,x 2 ). (5) 

The advantage for working with T± is that the non-singlet parts of T± does not mix under renormalization. 
The singlet parts are certainly mixed with the purely gluonic matrix elements. We only need to consider 
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the scale dependence of T + . The scale dependence of T_ can be obtained from the properties given in 
Eq.(|2|). At one-loop level, we can divide the scale-dependence or the evolution into three parts: 

dT + (xi,x 2 , n) _ dT + (x 1 ,x 2 ,fi) 
dlrifi dln.fi 

In the above the first two parts are of the non-singlet part. The first part is a convolution of T + (xi,X2) 
with xip > 0. The second part is a convolution with T + (xi, x 2 ) with x\ < or x 2 < 0. The third part is 
the mixing part with the purely gluonic matrix elements. Detailed expressions of each part will be given. 

Before we study the scale dependence from each part we briefly explain our method. Our method is 
to directly calculate these twist-3 matrix elements with parton states instead of a hadron state. If we use 
helicity instead of the spin vector s to describe the spin, the above defined matrix elements with twist-3 
operators, generically denoted as O, are the non-diagonal parts of the matrix elements {P, \'\0\P, X) in 
the 2x2 helcity space. That is we need to study the forward scattering amplitudes with helicity-flip. 

If we use a single quark to replace the hadron, we will always have T± = because helicity conservation 
of QCD. However, we can use a multi-parton state instead of a single-quark state. Factorizations of SSA 
have been studied in [I2j HU [TU [15] by using multi-parton states. We can introduce the following 
multi-parton state: 

\n[X}) = \q(p, X q )[X = Xg\) + C™\q( Pu X q )g(p 2 , X g )[X = X q + X g }) + • • • , (7) 

with pi + p 2 = p. The qg-st&te is in the fundamental representation of 5f7(iV c )-gauge group. We take 
the momenta as p\ = xop and p 2 = (1 — xq) P with p^ = P^ = (P + , 0, 0, 0). C q9 is a real constant. If we 
calculate T± with this state in Eq.([T]), we will find that T± receives nonzero contributions only from the 
matrix elements of the interference between the single quark- and the gg-state, i.e., (q(X q )\0\q(X q )g(X g )} 
or (q(X q )g(X g )\0\q(X q )) . It is noted that the total helicity in the bra- and ket-state is different, but the 
quark always has the same helicity. At tree-level, we obtain: 

r| 0) (xi,x 2 ) = g s TTC q9 V2x^(N^ -l)(x 2 -x 1 )5(l-x 1 )5(x 2 -x ), 

T {0) (x 1 ,x 2 ) = -gsirC^V^iN^ -l)(x 2 -x 1 )5(l-x 2 )5(x 1 -x ). (8) 

Because the multi-parton state depends on xq with xo < 1 per definition, T± calculated with the state 
will also depend on xq. But the scale dependence or the evolution will not depnd on xo and C qg . 

It is possible to have more multi-parton states in Eq.([7|), represented with However, not all 
possible multi-parton states for the cases with twist-3 operators are needed [H]. The basic idea for cases 
with twist-3 operators is to consider various helicity-flip matrix elements like {a,b\0\c) and (c\0\a,b) 
with a, b and c as possible partons in QCD. Because of helicity flip we only need to consider those 
matrix elements with three combinations of partons: {q,g\0\q), (q, q\0\g) and {g, g\0\g). The complex 
conjugated matrix elements should also be included for consistency. Because of the three combinations 
we can divide the evolution in Eq.([6]) into three parts at one-loop. For each combination with a given 
operator O one can construct the corresponding 2x2 spin density matrix in helcity space for a spin-1/2 
system. The non-diagonal part is relevant for the transverse polarization. For O being those operators 
used to define twist-3 matrix elements, one can extract these matrix elements from the corresponding 
non-diagonal parts of spin density matrices. The spin-density matrices for each combination are given in 
[H]. We will also call the contribution extracted from (q,g\0\q), (q, q\0\g) and (g, g\0\g) as qg-, qq- and 
gg contribution, respectively. For detailed description of those parton states and spin-density matrices 
we refer to |14j . We will also use the same notations for these multi-parton states as used in |14j . 
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dT + (x 1 ,x 2 ,fi) 
d\nu 
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dT + (xi,x 2 , fi) 
dlnu 



(6) 
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Figure 1: A set of diagrams of one-loop corrections to T±(x\,x 2 ) from the gg-contribution. This set only 
contains the self-energy corrections represented by black dots, and corrections with one gluon emission 
from a gauge link. 

Now we turn to the first part in Eq.([6]). The scale dependence of this part can be written as the 
evolution: 



d T±(x 1 ,x 2 ,n) 



</.</ 



a f 1 

— / d£id£2J r ±{x 1 ,X2,€i,€2)T±(£i,t2,fJ>). (9) 
7T Jo 



With the tree-level result given in Eq.([8|). the kernel J-+ can be determined by calculating the one- loop 
correction of T + . The one-loop contributions, denoted as T+\xi,X2,xq), are represented by diagrams 
given in Fig.l and Fig. 2. With the tree-level result and the scaling property of we have at the leading 
order: 

where • • • denote the contribution from the ^t-dependence of g s appearing in the definition in Eq.(l). We 
notice here that the calculation of a?2/£i> £2/^1) > hence of F + , is slightly different than that 

of r} 1} (xi ,x 2 ,x ). In the later, all variables £0,1,2 are always smaller than 1. Especially, with x < 1 one 
finds that T + (xi,x 2 ,Xo) are nonzero only for x\ t 2 < 1 and \x\ — x 2 \ < 1. While in the former, any of the 

variables in T+\xi/£\,X2/£i ,£2/^1) can t> e larger than 1. We take Fig.lh as an example to explain the 
difference. 

We denote the momentum of the gluon emitted from the gauge link as k, the contribution to 
(xi , X2 , Xq) from Fig.lh can be written in the form: 



T + (xi,x 2 ,x ) 



lh 



1 elk d k_\_ 1 1 r , , , 

X M ' - X ^ I (2vr)3 fe 2 + ie(p2 _ fe) 2 + fe [- ' ] (U) 
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with k + fixed as p^ — k + = (1 — X2)p + ■ We use the dimensional regularization with d = 4 — e. /i, is the 
renormalization scale. There are two poles from the two propagators in the /c~-plane. Their positions 
are determined by 1 — x<± and X2 — xq, respectively. Performing the k~- and /^-integration one has the 
nonzero contribution proportional to 5(1 — X\)9(x2 — xq)9(1 — X2) to T+(x\, X2, xq). From the above 
expression for one can also obtain the kernel before the integrations by taking the derivative against 
In /i. We have then 



lit 



oc lim 5(1 — xi/£i)e// / 
e->o J 



dk~d d - 2 k ± 



(2vr)3 2(s 2 /£i - Z 2 /ti)p + k- -k 2 ± + ie 



-2(l-x 2 /£i)p + k- -k\ + ie l 



:i2) 



Now we can see that the positions of the two poles are now determined by (1 — a?i/£i) and (x 2 /£i —^2/^1)-, 
respectively. Therefore, in the calculation of the kernel the positions of poles in the /c~-plane are differently 
determined than those in the calculation of T + (xi, X2, xq). It is then straightforward to obtain: 



Ih 



Xl 



Xl)[9(x2-i2)e(x 1 

x 2 2 2 



X 2 ) - 0(& - X2)6(X2 - Xl) 



(xi 



6) 2 



X2 - & 



+ 



Xi 



(13) 



We notice here that there is a singularity in the second term in [• • •] for £2 ~ xi. This is a light-cone 

singularity which will be canceled. We also note that the first term is nonzero as a distribution for 
x\ — > X2, although it is proportional to x\ — X2- 
The result of Fig. If reads: 



J r + (xi,x 2 ,£i,&, 



1/ 



<K£l - xi)6(& - x 2 ) 



Ki-xal dy 

y 



(14) 



This contribution also contains a light-cone singularity. This will be canceled by that in the contribution 
from Fig.lh. In the final result there is no light-cone singularity. In calculating the above contribution a 
special care should be taken. At first look the integration variable y may be changed by y — > y/\xi — x 2 \ 
so that one obtains the integral over y from to 1 instead of to \x\ — x 2 \. But this change is not 
correct. The correct result is only obtained by keeping the interpretation of y as a momentum fraction in 
unit of p + . This can be shown as in the following: The origin of the integral comes from the integration 
over k + with the eikonal propagator l/(k + + ie), i.e., the propagator from the gauge link represented 
by the double line in Fig. If. The light-cone singularity comes from the region with k + ~ 0. One can 
regularized this singularity in such a propagator with the so-call rj- or A-regulator used in |16t [T7] . With 
the regularization one can calculate all contributions involving the eikonal propagator. The results at 
end are the same as we have done for Fig. If by keeping y as the momentum fraction in unit of p + . 

We notice that the contribution from Fig. lb introduces an A^^-dependence in the kernel with Np as 
the number of quark flavors. By taking the scale dependence of g s in the definition in Eq.(pQ) into account, 
the iVi?-dependence is canceled. Adding all contributions from Fig.l to J-+ and the contribution from 
the scale dependence of g s in the definition we have: 



J r +(x l ,x 2 ,( l ,^ 2 ) 



Fig.l 



X2) 



-N c In \xi — x 2 \ + 



N 2 

c 



l-llnxa-ilnx! 
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+^<K6 " xi) ( " x 2 )0(x2 ~ 6) - 0(x 2 - xi)0(6 - x 2 ) 



Xl - x 2 



(xi - 6) 2 («2 - 6)+ si - 6 

+ ^<5(6 - X 2 ) (0( Xl - X 2 )0(6 - Xl) - 0(x 2 - X!)^(X1 - 6) 



Xl - x 2 



+ 



(6~X 2 ) 2 6~X 2 



' 2 



x 2 



<5(6-xi)0(6-x 2 )-± 



1 



1 



6V(6-x 2 ) + 6-xi 



+5(6-x 2 )0(6-xi)^ 



1 



1 



6 V(6 - xi)+ 6 - X2 



L. - x 2 - a + 6)^(6 - x 2 ) ( ? + ? 



2iV c 

The +-distributions appearing in the above are defined as: 



1 



v6 6/(6-x 2 ) 







(x - 0+ 







dq h / (x) mx, 



x — £ 

/(0 - f_M 

£ — X 



+ /(x)ln(l-x) 



These +-distributions are different than the standard +-distribution which will be given later. 
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Figure 2: Another set of diagrams for one-loop corrections of T± from the gg-contribution. 

The contributions from Fig. 2 can easily be worked out. The results are: 
N c 5(6 - x 2 ) 



J r +(xi,x 2 ,6,6^ 



j" + (xi,x 2 ,6,6) 



2d 



6-6 



0(6 - xi) 



2^^(x 2 -xi) 
6 



0(xi - X 2 ) f Xl 

— 7 t — 2— x 2 - x 2 - Xl 

6-6 V 6 



0(X1 - 6)0(^2 - Xl) 



Xl - x 2 

6 - x 2 



N, 



5(6 -xi) 



26 



fl(x 2 -xi)g(6~x 2 ) 

6(6 -xi) 2 



(2x 2 - 2xix 2 - 2£ 2 - 6x2 + 3xi6 



6 



6(x 1 -x 2 ) { 2x 2 . x 1 -x 2 / .( x 2 6 

+— t — t- (xi - x 2 - 6)0(6 - x 2 + t-0 (a? a - 6) 3 - 2— - 2— 

v«i&2 xi-6 V 6 6, 



2c 



2r/ 



2iV c 

£(6 - xi) x 2 - xi + g 2 
2A^c x x (6-xi) 
xi - x 2 



2/ 



6 

2N C 



-0(x 2 -xi +6) + 



£(6 - xi)0(xi - x 2 ) 



9{xi - x 2 ) 

x 2 



xi - 6 

x 2 (xi - x 2 ) 
x?(6 - XI ) ' 



66 ' 



x\ -x 2 - 6 



(17) 



the contributions from Fig.2e, 2g, 2h and 2i are zero. We do not try to give the sum these contributions 
because the result is too lengthy. Our final result for J 7 ^, i.e, the first part is then given by Eq.([9]) with 
the kernel: 

(18) 



F+=F+ + 

^ifl-l i=a,b,c,d,f 

In this part it is trivial that the kernel does not depend on C q9 . 
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Figure 3: One-loop diagrams in the light-cone gauge for T±(xi,x 2 ) from the (/(/-contribution. 

Now we turn to the second part or the (/(/-contribution in Eq.([6j). As discussed before, in this case 
we need to consider the matrix element as (q(pi), q(p 2 ) \ 0\g{p)), where momenta of partons are specified. 
In constructing the corresponding spin density matrices for this part, one notes that the (/(/-state must 
have the total helicity A = 0. There are two possible (/(/-states with A = 0. One is the symmetric state in 
helicity state. Another one is antisymmetric state. Therefore, the general qq -state is a superposition of 
these two states. We denote the weight of the symmetric state in the superposition as C+ 9 , and the weight 
of the anitysmmetric state as C 99 - The qq is in color-octet. The tree- level result for the (/(/-contribution 
reads: 



T[°\x u x 2 ) = ng s {N* - l)x/2^xo 



(Cf - C q J)5( Xl + x )<5(x 2 - x ) + (Cf + C q J)5(x 2 + x )5( Xl - x ) 



.(19) 



At tree-level, T + (x\,x 2 ) of the (/(/-contribution is zero for xi >2 > 0, but nonzero for x\ < or x 2 < 0. 

At one-loop, T + (xi,X2) with x% j2 > can be nonzero. It receives contributions from diagrams in 
Fig. 3. In Fig. 3 the diagrams are for contributions in the light-cone gauge n ■ G = 0. In Feynman gauge 
there are more diagrams. The scale dependence of this part can be written as the convolution: 



d r + (xi,x 2 , aO 



qq 



— [ d&d&r+x (xi,x 2 , 6,6)0(6 6,6) 

vr Jo 
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+— / di 1 d^F +2 { Xl ,x 2 ,i 1 ,i 2 )e^ 2 -&)t + (6,6-6)- 

In the above we have 6,2 > 0. Prom Fig. 3 we can obtain the kernels as: 

1 X\ 

—<5(6 - 2:2)0(6 - xi)-2, 

SI 



(20) 



7+1(2:1,2:2,6,6) 
7+2(2:1,2:2, £1,62) 



2N r 



N c x , c w , (x 2 - Xl )(^ - xj + Xl x 2 ) 

"-^(6 - aJijeCsj - X 2 ) -2—7- ; r . 

2 62^1(6 - 2:1 + x 2 ) 



+7^-5(6-2:1)0(6-2:2)11 

, N c- 1 s:^ snf ,X 2 (x 2 -Xi) 

H — f^rr— <H6 - xi)6»(xi - s 2 )- 



^ 2 -* 1 )%^ + 0(*i-* 2 ) X2fe - X2) 



6 - 2:1 



•'■1 



2N r 



6*1 



(21) 



The derived kernels do not depend on the introduced weights C± as expected. 

For the mixing part with the purely gluonic operators one needs to calculate the contribution T + 
from the matrix element like {g(pi), g(p 2 )\0\g(p)) . The state with the two gluons must have the same 
color of the state with one gluon. Therefore, the colors of the two gluons can be coupled through if abc 
and d abc . The state with the two gluons must have the total helicity A = for helicity flip. Again, there 
are two possible states. The general two gluon state with A = is a superposition of the states. One is 
the symmetric state in helicity space, for which we introduce a weight J 7 ^ 9 for the color structure with 
if abc . Another is the antisymmetric state in helicity space, for which we introduce a weight T 9 ^ for the 
color structure with if abc . One can introduce weights in a similar way for the color structure with d abc . 
The corresponding spin density matrices for our calculation are defined in 




Figure 4: The one-loop diagrams for T±(xi,x 2 ) from the gg-contribution. The gluons from the bottom 
are labeled as, the left gluon is with pi, the middle one is with p 2 , both are incoming. The right one is 
with p as an outgoing gluon. 

(f) if) 

At tree-level, we obtain the results for T G ' and T^ G with x\ >2 > 0: 



T G f '°\x 1 ,x 2 ) = & c (iY c 2 -l) Xl x 2 (x 1 -x 2 ) 



V2 



C+ U(l - x 1 )5(x 2 - x ) - 5(1 - x 2 )5(xi - xo)^j 



-C- \5(1 - x\)8(x 2 - x ) - 5(1 - x 2 )5(xi - x ) 



T 



(/,o) 



V2 



A^ c (A" c 2 -l)xix 2 (xi-x 2 ) 



C + (s(l - xi)5(x 2 - x ) + 5(1 - x 2 )5(x\ - x ) 



-C_ (s(l - xi)5(x 2 - x ) + 5(1 - x 2 )5(xi - x )^j 
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c 4 



■99 



(22) 



For other cases, i.e., x\ %2 < the functions can also be obtained. Similar results for the color structure 
with d abc can be obtained with the replacement of color factors. At tree-level, T + is zero here, but nonzero 
at one-loop. At one-loop there are six diagrams which give contributions to X+. Two of them are given in 
Fig. 4. Other four diagrams can be obtained through permutations. Only the diagrams given in Fig. 4 will 
contribute to the mixing part in the scale dependence for x\p > 0. The calculations are straightforward. 
We find that the contributions from the two diagrams are renormailized with the purely gluonic matrix 
elements with x\ 2 > 0. To give our result we introduce the notation: 



T± G (x 1 ,x 2 ) = (T^\x 1 ,x 2 ) + T^(x 1 ,x 2 )) ± {t { / ) g (x 1 ,x 2 ) +T^ G ( Xl ,x 2 , 

The mixing- or third part of the scale dependence in Eq.(|6]) can be given as: 

dT+(xi,x 2 ,n) a s f 1 At * , c c \r u c\ 

- — / d£id£ 2 J r G {xi,X2,£i,&)T +G {€i,&), 



(23) 



(24) 



mi 



with 



FG(xi,X2,£l,& 



X\%2 



-b{x x -X2-C1 + &)0 (£1 ~ xi)—^ 

Z S2S1 



-5(x 2 -xi-£i + £2)9 (€1 ~ X2 



(6 



X2 



2&? 



(25) 



Again, the kernel derived here does not depend on the states, i.e., on the weights like J-±. 

With the three parts in Eq.Q given in the above, we have derived the scale dependence at one-loop 
for the twist-3 matrix element T + (x\, x 2 , /•*) for x% j2 > 0. From the derived scale dependence one can 
easily find the scale dependence of T_(x\, x 2 , fJ,) by using T_(xi, x 2 , ft) = T + (x 2 ,Xi, (j,). Two special 
cases can be derived from our general result. One is the evolution of the so-called soft-gluon-pole matrix 
element, which is obtained by taking x% = x 2 . We can derive from Eq.([6|) the evolution of T + (x,x, fi) by 
taking the limits x% — > X2 ± + . The results are same in these limits. We have with z = x/£: 



8T F (x, x, /i) 



OLs 
7T 



1 dz 

z 



P qg (z)T F (£,0 + 



N c (l + z)T F (x^)-(l + z 2 )T F (tO 



2N r 



2 1 

(1 - 2z)T F (x, x — £) + Ta jF (x, x — £) 



+ T a ,f(x,£) 



N c T F (x,x) 



dz(l- z) 2 + Z 2 



where the quark splitting kernel is 



Nt 



2N r 



1 + z 2 



and the standard +- distribution is defined as: 



dz- 



1 /(*)-/(!) 



/(l)ln(l 



(26) 



(27) 



(28) 
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Since we have derived the same result in Eq,(|26p from the two limits x\ — > X2±0 + , it indicates that there 
is no scale dependence of T& : p(x,x). Comparing existing results in (6j[7], we find that our result of the 
non-singlet part, i.e., the contributions in the first- and second line in Eq. ([26|) . is in agreement with that 
in [7]. The main difference between the results in [BJITIIEIE] is the last term in the second line. We notice 
that during preparing our work the last term is also confirmed in [18] . The mixing part is in agreement 
with that given in [6] by taking the difference of definitions into account. Our method for deriving the 
scale dependence is different than that in [U El El EJ HE]. Therefore, our result for the soft-gluon-pole 
matrix element gives an independent verification of existing results. 

Another special case is the evolution of soft-quark-pole matrix elements. In this case one has either 
x\ = or X2 = 0. Form our general results we can derive: 



dT + (0,x,fi) 



a s ( f 1 dz 
tt \Jx z 

3(n? - n , , f 1 dz 

+ v ; > t + o,x + / - 

4N r J x z 



2N C (l-z) + 2 (1-z). 



N 7 2 1 



dT + (x,0,n) 



a, 

7T 
+ 



1 1 f 1 dz tJ _ ^ A 

1 T + ((,i-x 



dz 



2N C (l-z)+ 2N C +V 2 (l-z)+ 



3(iV c 2 



AN r 



I . , f 1 dz 
-T + (x,0 + / — 

Jx Z 

l r 1 dz 



2 (1 - z) 



A ' ' T+(x,x-£)--Ut + (-£,0) 



2N r 



(29) 



It is interesting to note that the first- and second part in Eq.([6]) are not finite for x\ = or x<i = 0. They 
contain light-cone singularities. But they are canceled in the sum. 

To summarize: We have derived the evolution of the matrix elements of twist-3 quark-gluon operators. 
These matrix elements are important ingredients for predict SSA in QCD collinear factorization. From 
our results we have derived evolutions in two special cases. One is of soft-gluon-pole matrix elements, 
another is of soft-quark-pole matrix elements. A comparison has been made with the existing results 
for soft-gluon-pole matrix elements. In this letter, we have only considered the case with x\ 2 > in 
the evolution. The study of other cases including the evolution of purely gluonic matrix elements are in 
progress. 
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